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Abstract 

A study of energy behavior of the pion spectra and interferometry scales is carried out for the top SPS, 
RHIC and LHC energies within the hydrokinetic approach. The latter allows one to describe evolution of 
quark-gluon and hadron matter as well as continuous particle emission from the fluid in agreement with the 
underlying kinetic equations. The main mechanisms that lead to the paradoxical, at first sight, behavior of 
the interferometry scales, are exposed. In particular, a slow decrease and apparent saturation of R ou t / ' Rside 
ratio with an energy growth happens due to a strengthening of positive correlations between space and time 
positions of pions emitted at the radial periphery of the system. Such en effect is a consequence of the 
two factors: a developing of the pre-thermal collective transverse flows and an increase of the initial energy 
density in the fireball. 
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1. Introduction 

The pion interferometry method, which is based 
on the Bose-Einstein or Fermi-Dirac interference of 
identical particles, has been proposed first in Refs. 
[H I3] for measurements of the geometrical sizes 
and shapes of the interaction region in hadronic 
collisions. Then it has been developed in Refs. 
[S HI S HI as a tool for a study of rapidly expanding 
fireballs formed in ultra-relativistic heavy ion col- 
lisions. Despite the extremely small sizes of such 
systems - the order of the value is around 10 -14 
m - they have, nevertheless, a pronounced inhomo- 
geneous structure. The generalized treatment of 
the interferometry measurements asserts that the 
measured scales - the interferometry radii - are as- 
sociated iust with the homogeneity lengths in the 
system 7]. Only in the very particular case of a 
finite homogeneous system such lengths correspond 
to the total geometrical sizes, but normally they are 
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smaller than the latter. The interferometry scan- 
ning of femto-systems at various total momenta of 
pion pairs allows one to analyze the homogeneity 
lengths related to different space-time regions of the 
expanding fireball 0, Q . An understanding of in- 
terferometry in terms of the homogeneity length, 
as opposed to the simple-mind geometrical picture, 
provides explanation to some, at first sight para- 
doxical, results at RHIC. 

Naively it was expected that when the energy of 
colliding nuclei increases, the pion interferometry 
volume Vi n t - product of the interferometry radii 
in three orthogonal directions - rises at the same 
maximal centrality for Pb+Pb and Au+Au colli- 
sions just proportionally to However, when 
experiments at RHIC starts, an increase of the in- 
terferometry volume with energy turns out to be es- 
sentially flatter then the proportionality low. This 
is one of the components of so-called RHIC HBT 
puzzle 

Another item of the HBT puzzle is an unexpected 
RHIC result as for the ratio of the two transverse 
interferometry radii: one that is measured in the 
direction of the sum of transverse momenta of the 
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(pion) pair, R ou t, and another, R S ide, - orthogo- 
nally to it and the beam axis. While the latter is 
associated with transverse homogeneity length, the 
former includes besides that also additional contri- 
butions, in particular, the one which is related to a 
duration of the pion emission. Since the lifetime of 
the systems obviously should grow with collision en- 
ergy, if it is accompanied by an increase of the initial 
energy density and/or by a softening of the equa- 
tion of state due to phase transition between hadron 
matter and quark-gluon plasma (QGP), then the 
duration of pion emission should grow with en- 
ergy and Rout/Rside ratio could increase. A very 
large value of the ratio was predicted as a signal of 
the QGP formation The RHIC experiments 

brought the result: the ratio Rout/Rside ~ 1 and 
similar or even smaller than at SPS. 

In this letter we analyze the pion spectra and the 
intcrferometry scales as depending on the energy 
of central Pb+Pb and Au+Au collisions supposing 
the formation of QGP. The aims of the study are to 
describe these pion observables at the top SPS and 
RHIC energies and make predictions for LHC en- 
ergies, also to understand the physical mechanisms 
responsible for the peculiarities of energy depen- 
dence of the interferometry radii, in particular, the 
Rout/Rside ratio. Numerical calculations are based 



on the HydroKinetic Model (HKM) [ll|, |l2j that 
allows one to describe the evolution of quark-gluon 
and hadron matter and continuous particle emission 
from the fluid in agreement with the underlying ki- 
netic equations. 

2. Hydro-kinetic approach to A+A collisions 

Let us briefly describe the main features of the 
HKM [l2| ■ It incorporates hydrodynamical ex- 
pansion of the systems formed in A + A collisions 
and their dynamical decoupling described by escape 
probabilities. The method corresponds to a gener- 
alized relaxation time approximation for the par- 
ticle emission functions applied to inhomogeneous 
systems expanding into vacuum according to the 
Boltzmann equations. This method also allows one 
to describe the viscous effects at the hadronic stage 
of the evolution. The basic hydro-kinetic code, pro- 
posed in [l2| . is modified now to include decays of 
resonances into expanding hadronic non-chemically 
equilibrated system and, based on the resulting 
composition of the hadron-resonance gas at each 
space-time point, to calculate the equation of state 
(EoS) in a vicinity of this point. The obtained local 



EoS allows one to determine the further evolution 
of the considered fluid elements. The complete pic- 
ture of the physical processes in central Pb + Pb 
and Au + Au collisions encoded in calculations is 
the following. 

2.1. Initial conditions 

Our results are all related to the central rapidity 
slice where we use the boost-invariant Bjorken-like 
initial condition. We consider the proper time of 
thermalization of quark-gluon matter as the mini- 
mal one discussed in the literature, tq = 1 fm/c. 
The initial energy density in the transverse plane 
is supposed to be Glauber-like [l3j], i.e. is propor- 
tional to the participant nucleon density for Pb+Pb 
(SPS) and Au+Au (RHIC, LHC) collisions with 
zero impact parameter. The height of the distri- 
bution - the maximal initial energy density - e(r = 
0) = eo is the fitting parameter. From analysis of 
pion transverse spectra we choose it for the top SPS 
energy to be e = 9 GeV/fm 3 ((e} = 6.4 GeV/fm 3 ), 
for the top RHIC energy e = 16.5 GeV/fm 3 ((e) 
= 11.6 GeV/fm 3 ). The brackets < ... > correspond 
to mean value over the distribution associated with 
the Glauber transverse profile. We also demon- 
strate results at eo = 40 GeV/fm 3 and eq — 60 
GeV/fm 3 that probably can correspond to the LHC 
energies of 4 and 5.6 ATeV. As for the latter, ac- 
cording to at the top LHC energy e max re- 
calculated to the time 1 fm/c is 0.07 * 700 = 49 
GeV/fm 3 . We suppose that soon after thermaliza- 
tion the matter created in A+A collision at RHIC 
energies is in the quark-gluon plasma (QGP) state. 
Also, at time Ti, there is peripheral region with rel- 
atively small initial energy densities: e(r) < 0.5 
GeV/fm 3 . This part of the matter ("corona") does 
not transform into QGP and have no chance to be 
involved in thermalization process (ljj . By itself 
the corona gives no essential contribution to the 
hadron spectra [l5| . This part of the matter, prob- 
ably, is not thermalized. One should consider it 
separately from the thermal bulk of the matter and 
should not include in hydrodynamic evolution. 

At the time of thermalization, To = 1 fm/c, the 
system already has developed collective transverse 
velocities 

nana. 

The initial transverse rapidity 
profile is supposed to be linear in radius rx- 



ijt = a -pr~ where Rt 
Rt 



< r T >., 



(1) 



Here a is the second fitting parameter. Note that 
the fitting parameter a should include also a pos- 
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itive correction for underestimated resulting trans- 
verse flow since in this work we do not take di- 
rectly into account the viscosity effects [l8| neither 
at QGP stage nor at hadronic one. In formalism 
of HKM 12| the viscosity effects at hadronic stage 
are incorporated in the mechanisms of the back re- 
action of particle emission on hydrodynamic evolu- 
tion which we ignore in current calculations. Since 
the corrections to transverse flows which depend on 
unknown viscosity coefficients are unknown, we use 
fitting parameter a to describe the " additional un- 
known portion" of flows, caused both factors: by a 
developing of the pre-thermal flows and the viscos- 
ity effects in quark-gluon plasma. The best fits of 
the pion transverse spectra at SPS and RHIC are 
provided at a = 0.194 {{v T ) = 0.178) for SPS en- 
ergies and a = 0.28 ((v T ) = 0.25) for RHIC ones. 
The latter value we use also for LHC energies aim- 
ing to analyze just influence of energy density in- 
crease. From a fitting of the pion spectra we can 
conclude that an " additional portion" of flow is big- 
ger at RHIC than at SPS where the parameter a 
is com pat ible with that free streaming of partons 
brings jig, [l7|. One can conclude, thus, that longer 
viscous QGP evolution at RHIC has a bigger im- 
pact for a developing of transverse flows than at 
the SPS energies. 

2.2. Equation of state 

Following Ref. [l2| we use at high tempera- 
tures the EoS [l|| adjusted to the QCD lattice 
data with the baryonic chemical potential fiB = 
and matched with chemically equilibrated multi- 
component hadron resonance gas at T — 175 MeV. 
Such an EoS could be a good approximation for the 
RHIC and LHC energies; as for the SPS energies we 
utilize it just to demonstrate the energy- dependent 
mechanism of formation of the space-time scales. 
We suppose the chemical freeze-out for the hadron 
gas at T c h — 165 MeV [20]. It guarantees us the 
correct particle number ratios for all quasi-stable 
particles (here we calculate only pion observables) 
at least for RHIC. Below T c h a composition of the 
hadron gas is changed only due to resonance de- 
cays into expanding fluid. We include 359 hadron 
states made of u, d, s quarks with masses up to 2.6 
GeV. The EoS in this non-chemically equilibrated 
system depends now on particle number densities 
rii of all the 359 particle species i: p — p(e, {rii}). 
Since the energy densities in expanding system do 
not directly correlate with resonance decays, all the 
variables in the EoS depend on space-time points 



and so an evaluation of the EoS is incorporated in 
the hydrodynamic code. We calculate the EoS be- 
low T c h in the Boltzmann approximation of ideal 
multi-component hadron gas. 

2.3. Evolution 

At the temperatures higher than T c h the hydro- 
dynamic evolution is related to the quark-gluon and 
hadron phases which are in chemical equilibrium 
with zero baryonic chemical potential. The evolu- 
tion is described by the conservation law for the 
energy-momentum tensor of perfect fluid: 



d v T> iV {x) = 



(2) 



At T < T c /,.=165 MeV the system evolves as non- 
chemically equilibrated hadronic gas. The con- 
cept of the chemical freeze-out implies that after- 
wards only elastic collisions and resonance decays 
take place because of relatively small densities al- 
lied with a fast rate of expansion at the last stage. 
Thus, in addition to (0), the equations accounting 
the particle number conservation and resonance de- 
cays are added. If one neglects the thermal motion 
of heavy resonances the equations for particle den- 
sities ni(x) take the form: 

d fl (rii(x)u fJ '(x)) = -Tini(x)+^2bijTjnj(x)(3) 

3 

where by — BijNij denote the average number of 
ith particles coming from arbitrary decay of jth res- 
onance, Bij — Tij/Tj t tot is branching ratio, Nij is 
a number of ith particles produced in j — ^ i decay 
channel. We also can account for recombination in 
the processes of resonance decays into expanding 
medium just by utilizing the effective decay width 
Fi, e ff — 7IV We use 7 = 0.75 supposing thus that 
near 30% of resonances are recombining during the 
evolution. All equations ([2]) and 359 equations ([3]) 
are solved simultaneously with calculation of the 
EoS, p(x) = p(e(x), {ni(x)}), at each point x. 

2.4- System's decoupling and spectra formation 

During the matter evolution, in fact, at T < 
T c h, hadrons continuously leave the system. Such 
a process is described by means of the emis- 
sion function S(x,p) which is expressed for pions 
through the gain term, G 7r {x,p) 1 in Boltzmann 

equations and the escape probabilities V- K {x,p) = 

00 

exp(- / dsRir+his^r + ^(s - t),p)): S- !T {x,p) = 
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G 7r (x,p)V 7r (x^p) [h, 12 1 . For pion emission in re- 
laxation time approximation G-ir ~ firRir+h + GH^Tr 
where f n (x,p) is the pion Bose- Einstein phase- 
space distribution, R 7r+ h(x,p) is the total collision 
rate of the pion, carrying momentum p, with all 
the hadrons h in the system in a vicinity of point 
x, the term Gh^-k describes an inflow of the pions 
into phase-space point (x,p) due to the resonance 
decays. It is calculated according to the kinematics 
of decays with simplification that the spectral func- 
tion of the resonance H is S(p 2 — (mi/) 2 ). The cross- 
sections in the hadronic gas, that determine via the 
collision rate R^+h the escape probabilities V(x,p) 
and emission function S(x,p), are calculated in ac- 
cordance with the UrQMD method [2l[ . The spec- 
tra and correlation functions are found from the 
emission function S in the standard way (see, e.g., 



11]). 



3. Results and conclusions 

The pion emission function per unit (central) ra- 
pidity, integrated over azimuthal angular and trans- 
verse momenta, is presented in Fig. 1 for the top 
SPS, RHIC and LHC energies as a function of trans- 
verse radius r and proper time r. The two fitting 
parameters eo and (vt) are fixed as discussed above 
and marked in figures. The pion transverse mo- 
mentum spectrum, its slope as well as the absolute 
value, and the interferometry radii, including R ou t 
to R s ide ratio, are in a good agreement with the 
experimental data both for the top SPS and RHIC 
energies. 

As one can see particle emission lasts a total life- 
time of the fireballs; in the central part, r « 0, 
the duration is half of the lifetime. Neverthe- 
less, according to very recent results [l2|, 22 1, the 



Landau/Cooper- Frye presentation of sudden freeze- 
out could be applied in a generalized form account- 
ing for momentum dependence of the freeze-out 
hypersurface cr p (x); now o~ p {x) corresponds to the 
maximum of emission function S{t tT {v,p) 1 T^p) at 
fixed momentum p in an appropriate region of r. 
This finding gives one possibility to keep in mind 
the known results based on the Cooper-Frye for- 
malism, applying them to a surface of the maximal 
emission for given p. Then the typical features of 
the energy dependence can be understood as fol- 
lows. The inverse of the spectra slopes, T e //, grows 
with energy, since as one sees from the emission 
functions, the duration of expansion increases with 



initial energy density and, therefore, the fluid ele- 
ment gets more transverse collective velocities vt 
when reach a decoupling energy densities. There- 
fore the blue shift of the spectra becomes stronger. 
A rise of the transverse collective flow with energy 
leads to some compensation of an increase of R S ide- 
qualitatively the homogeneity length at decoupling 
stage is R side = Rceom/y/^ + (v^)m T /2T, (see, 
e.g., 0). So, despite the significant increase of the 
transverse emission region, Rceom, seen in Fig.l, 
a magnification of collective flow partially compen- 
sates this. It leads to only a moderate increase of 
the R s ide with energy. Since the temperatures in 
the regions of the maximal emission decrease very 
slowly when initial energy density grows (e.g., the 
temperatures for SPS, RHIC and LHC are corre- 
spondingly 0.105, 0.103 a nd 0.95 MeV for p T = 0.3 
GeV/c ) the Ri ong ~ t^T/tut [j| grows propor- 
tionally to an increase of the proper time associated 
with the hypersurface o~ PT (x) of maximal emission. 
As we see from Fig. 1 this time grows quite mod- 
erate with the collision energy. 

A non-trivial result concerns the energy behav- 
ior of the Rout I Rside ratio. It slowly drops when 
energy grows and apparently is saturated at fairly 
high energies at the value close to unity (Fig.l). 
To clarify the physical reason of it let us make a 
simple half-quantitative analysis. As one can see in 
Fig. 1, the hypersurface of the maximal emission 
can be approximated as consisting of two parts: the 
"volume" emission (V) at r w const and "surface" 
emission (S). A similar picture within the Cooper- 
Frye prescription, which generalizes the blast-wave 
model [23| by means of including of the surface 
emission has been considered in Ref. [24|]. If the 
hypersurface of maximal emission f(r) is double- 
valued function, as in our case, then at some trans- 
verse momentum px the transverse spectra and 
HBT radii will be formed mostly by the two con- 
tributions from the different regions with the ho- 
mogeneity lengths \ t v = \J < (A?'i) 2 > (i = side, 
out) at the ^-hypersurface and with the homo- 
geneity lengths Xi t s at the S-hypersurface. Sim- 
ilar to Ref.[8|, one can apply at tht/T 3> 1 the 
saddle point method when calculate the single and 
two particle spectra using the boost-invariant mea- 
sures nv = dcr^pt 1 = T(r)rdrdcf)dri(mT cosh(?y — 



y) 



Pr^p- cos((/> - a)) and ^ s = da^p^ = 

r(T)TdTd<pdri(—mT cosh(7y — y) dr ^ +px cos(0 — a)) 
for V- and S- parts of freeze-out hypersurface cor- 
respondingly (here rj and y are space-time and par- 
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Figure 1: The p^-integrated emission functions of negative pions for the top SPS, RHIC and LHC energies (top); the inter- 
ferometry radii (middle) Rout/ Rside. ratio and transverse momentum spectra (bottom) of negative pions at different energy 
densities, all calculated in HKM model. The ex peri mental data are taken from CERES [25ll and NA-49 Collaborations |2(1 l27l| 
(SPS CERN), STAR [HLH and PHENIX [U \EM Collaborations (RHIC BNL) 



ticle pair rapidities, the similar correspondence is 



for ang les <fi an d ot, also note that 



> 



rfr(r) 



12t 122]). Then one can write, ignoring for sim- 
plicity the interference (cross-terms) between the 
surface and volume contributions, 

Rside = c V^side,V + c S^lide,S (4) 
Rout = c V^out,V + c S^OMt,s(l — d") 2 ' ^ 

where the coefficients c v + c| < 1 and we take 
into account that at p° /T 3> 1 for pions /3 ou t = 
Pout/p° ~ 1- All homogeneity lengths depend on 
mean transverse momentum of the pion pairs pt- 
The slope in the region of homogeneity ex- 
presses the strength of r— t correlations between the 
space and time points of particle emission at the S*- 
hypersurface f (r). The picture of emission in Fig. 1 



shows that when the energy grows the correlations 
between the time and radial points of the emission 
become positive, $f- > 0, and they increase with 
energy density. The positivity is caused by the ini- 
tial radial flows 16| u r (To), which are developed at 
the pre-thermal stage, and the strengthening of the 
t — t correlations happens because the non-central 
zth fluid elements, which produce after their expan- 
sion the surface emission, need more time Ti(eo) to 
reach the decoupling density if they initially have 
higher energy density eo- (Let us characterize this 

dTi(eg) 



effect by the parameter k 



den 



> 0). Then 



the fluid elements before their decays run up to 
larger radial freeze-out position rf. if a is the aver- 
age Lorentz-invariant acceleration of those fluid el- 
ements during the system expansion, then roughly 
for ith fluid elements which decays at time ti we 
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have at ar, > 1: r^n) « r i (r )+r 4 + (ii[(ro)-l)/a. 
Then the level of r — r correlations within the 
homogeneous freeze-out "surface" region, which is 
formed by the expanding matter that initially at 
To occupies the region between the transversal radii 
7"! (t ) and r 2 (r ) > ri(r ), is 

dr _ ri(ri) -r 2 (r 2 ) _ ^ _ R_ 
dr t\ - T 2 CqK 

and, therefore, the strength of r — r correlations 
grows with energy: 4- — >• 1. Note that here we 
account for r 2 — n « «(eo(r2(ro)) - eo(^i(ro))) and 
that de ^ r ^ ss — ^ where e = eo(r = 0) and i? is 
radius of nuclear. As a result the second S-term 
in Eq. ([5]) tends to zero at large eo , reducing, 
therefore, the R ou t/Rside ratio. In particular, if 
^"side v ^side s then, accounting for a similarity 
of the volume emission in our approximation and 
in the blast wave model, where as known \ S ide,v ~ 
A ut v, one can get: w 1 + consi ■ — s- 1 

at eo — > oo. It is worth noting that measure /xs 
also tends to zero when jj^ — > 1 that again reduces 
the surface contribution to side— and out— radii at 
large p T - 

The presented qualitative, in fact, analysis 
demonstrates the main mechanisms leading to the 
non-trivial behavior of R ou t to R s %de ratio exposed 
in detailed HKM calculations, see Fig.l (bottom, 
right). 

4. Summary 

We conclude that the energy behavior of the pion 
interferometry scales and transverse spectra can be 
understood if they are analyzed within fairly de- 
veloped hydro-kinetic models. The latter should 
be based on EoS which accounts for a crossover 
transition between quark-gluon and hadron mat- 
ters at high collision energies and non-chemically 
equilibrated expansion of the hadron-resonance gas 
at the later stage. The process of particle emis- 
sion from expanding fireball, that is not sudden and 
lasts about system's lifetime, should be correctly 
treated. The pre-thermal formation of transverse 
flow have to be taken into account. Then the main 
mechanisms that lead to the paradoxical behavior 
of the interferometry scales find a natural explana- 
tion. In particular, a slow decrease and apparent 
saturation of R ou t / Rside ratio around unity at high 
energy happens due to a strengthening of positive 
correlations between space and time positions of pi- 
ons emitted at the radial periphery of the system. 



Such an effect is a consequence of the two factors ac- 
companying an increase of collision energy: a devel- 
oping of the pre-thermal collective transverse flows 
and an increase of initial energy density in the fire- 
ball. 
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